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We apply a fundamental definition of time delay, as the difference between the time a particle 
spends within a finite region of a potential and the time a free particle spends in the same region, 
to determine results for photoionization of an electron by an extreme ultraviolet (XUV) laser field 
using numerical simulations on a grid. Our numerical results are in good agreement with those of 
the Wigner-Smith time delay, obtained as the derivative of the phase shift of the scattering wave 
packet with respect to its energy, for the short-range Yukawa potential. In case of the Coulomb 
potential we obtain time delays for any finite region, while - as expected - the results do not converge 
as the size of the region increases towards infinity. The impact of an ultrashort near-infrared probe 
pulse on the time delay is analyzed for both the Yukawa as well as the Coulomb potential and is 
found to be small for intensities below 10 13 W/cm 2 . 

PACS numbers: 33.80.Rv, 33.80.Wz 



I. INTRODUCTION 

The development of attosecond extreme ultraviolet 
(XUV) laser technology in recent years has offered the 
opportunity to observe and control the dynamics of elec- 
trons and the coupling to nuclear dynamics in atoms and 
molecules on their natural time scale. In particular, the 
capability to lock XUV pulses to a near-infrared (near- 
IR) pulse has initiated the development of techniques in 
which the dynamics is triggered by the attosecond pulse 
and observed as a function of the delay between the XUV 
and the near-IR pulses. Experimental observations in- 
clude, among others, the time resolution of the Auger de- 
cay [l| , the dynamics of electrons in valence shell [2| and 
excited states [H, H| , shake-up processes [H , and delays in 
the photoemission of electrons from different bands in a 
solid Q or different sub-shells in an atom [H, 0] ■ 

In particular, observations of substantial time delays 
during photoionization of atoms have generated signif- 
icant theoretical interest (e.g., [(1-1221]). These mea- 
surements are often analyzed in terms of the so-called 
Wigner-Smith (WS) time delay (e.g., 0-11, M, El GJ, 
2 Till). The WS time delay originally accounted for 
the delay between the propagation of a particle in a po- 
tential and that of a corresponding free particle over a 
distance R in an atomic or molecular scattering scenario 
23 . |2~H . In the latter context it has been pointed out 
24| , that this definition leads to a well-defined time de- 
lay as R — > oo as long as the potential vanishes quickly 
enough at large distances. Thus, for long-range poten- 
tials, such as the Coulomb potential, the asymptotic WS 
time delay is an intrinsically ill-defined concept. Hence, 
the related convenient calculation via the derivative of 
the phase shift of the wavefunction with respect to the 
energy of the particle appears to be not always applicable 
to the photoionization process. 

Having this in mind we note that an important aspect 
of the recent observations of time delays using attosecond 
XUV pump and near-IR probe pulses is the impact of 
the probing pulse on the result of the measurement 0- 



EIL 123- Since the probe pulse is usually ultrashort, the 
effect occurs over a finite time and, hence, during the 
propagation of the electron over a finite distance R in the 
potential only. Since the above mentioned time delay is 
a well-defined quantity over any finite distance, it can be 
therefore used to analyze the impact of the probing pulse 
for short- as well as long-range potentials. Furthermore, 
it is possible to study the effect using time-dependent 
numerical simulations on a grid, which intrinsically need 
to incorporate finite times and distances as well. 

In this article we report about such an analysis for 
photoionization of an electron in the short-range Yukawa 
potential as well as in the long-range Coulomb potential. 
To this end, we provide the basic definitions for the calcu- 
lation of a time delay with and without a strong probing 
field and discuss the application in numerical simulations. 
We then show that the numerical results for single XUV 
photoionization are well-defined over finite distances but 
clearly do not converge as R — > oo in the Coulomb case, 
in agreement with the discussion given in the early work 
by Smith [24[ . On the other hand, there is a good agree- 
ment between our numerical results and those for the 
original WS time delay in case of a short-range poten- 
tial. Finally, we investigate the impact of an ultrashort 
near-IR probing pulse on the time delay. Our results in- 
dicate that the effect is small as long as the intensity of 
the probe pulse does not exceed 10 13 W/cm 2 . 



II. THEORETICAL METHOD 



In this section we introduce the theoretical method, 
which we use to obtain time delays associated with the 
ionization of a target system in numerical simulations. 
To this end, we will first provide a set of basic definitions 
used in the method before we discuss its use in numerical 
simulations. 
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Basic definitions 



For a particle in a given state ^(r,t) the time spent in- 
side a region R with a potential V(r) can be expressed as 
(Hartrec atomic units, e = m = h = 1 are used through- 
out the paper): 



dt / dr\^(r,t)\' 



(1) 



While ty t R is, in general, finite for finite regions and any 
^(r, t), it is useful to compare f# ^ to the time spent by 
the free particle in R (or another reference time): 



t 



dt / dr|* (0) (r,i)p 



(2) 



as the "radius" of R grows to infinity [H| . Here, * (0) (r , t) 
is the free particle state corresponding to ty(r,t). The 
difference between tqi^R and i^(oj # defines the time delay 
associated with ^(r^i), the region R and the potential 
V{r): 



Ai$ o — r — i 



(3) 



The quantity Ai^j? is known to have a finite limit as 
the radius of R grows to infinity if the interaction van- 
ishes quickly enough [24| . Thus, At-qi ,r—>oo (and the as- 
sociated quantum mechanical operator) is well-defined 
for short-range potentials V(r) only. Provided that the 
limit exists, At^M^oo can be also expressed as the energy 
derivative of the phase shift tp induced by the potential 
V(r): 



dip 
dE 1 



(4) 



which is commonly known as the Wigner-Smith (WS) 
time delay. 

We intend to obtain time delays in the form of the 
time difference, given in Eq. ([3]), for a photoionization 
process. Physically, we are interested in the time that an 
- initially bound - electron needs to leave a certain re- 
gion (centered about the location of the residual target) 
following ionization due to the interaction with an exter- 
nal light field. To this end, we note that the expressions 
above can be readily applied to a particle in a superposi- 
tion of states and therefore consider the ionizing part (i.e. 

the continuum parts) of the wavefunction ^[ wn \r,t) in 
our adoption of Eq. ((T|) : 



J — oo J E 



dr — 

r C 
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(ion) 



dr\*[ ion) (r,t ^oo)| 2 



(5) 



We renormalize the wavefunction in order to be able to 
compare times and time delays arising for the ionization 
from different initial bound states ^i(r,t = 0). We can 
then define the time delay associated with the ionization 



from a specific initial bound state analogous to Eq. (|3]) 
as: 



At 



(6) 



where ^f\r,t) is the free particle state corresponding 
to the ionizing part of the wavefunction after transition 
from the initial state ^i(r, t = 0). According to this defi- 
nition we expect negative values for the time delays, since 
a free wave packet should spend more time in a given re- 
gion R than the corresponding wave packet that has the 
same asymptotic energy propagating in an attractive po- 
tential. We expect that Atq, it R has a well-defined finite 
limit (for R —> oo), i.e. the Wigner-Smith time delay, for 
short-range potentials, but not necessarily for long-range 
potentials such as the Coulomb interaction. In view of 
the intrinsical negative time delays for finite regions, we 
expect that the limit value is negative as well. We also 
consider the difference in the time delays for the ioniza- 
tions from two different initial states ^i(r, t — 0) and 
^j{r,t = 0) as: 



AT(*. t , R) = At* i:R - Aty.R. 



B. Numerical simulations of time delays 



(7) 



In order to use the above definitions in a numerical 
simulation of a photoionization process we need to iden- 
tify the ionizing part of the wavefunction ^f on '(r, t) as 

well as the corresponding free particle state ^l°\r,t). 
Since it is not straightforward to obtain e.g. the time of 
ionization (e.g. [§j]) and the form of the wave packet af- 
ter the transition into the continuum, it appears to be 
difficult to make use of Eq. ([5]) in a numerical simula- 
tion directly. We circumvent this obstacle by using the 
back-propagation technique. 

We first solve the time-dependent Schrodinger equa- 
tion (TDSE) of the system, initially in the state ^(r, t = 
0), under the interaction with the external light field on 
a space-time grid: 

i JU(r, *) = ( y + V(r) + (t) ) *(r, t), (8) 

where p is the momentum operator and Might (t) repre- 
sents the interaction with the ionizing light field. After 
the end of the interaction with the light field we separate 
the ionizing part of the wavefunction from the remaining 
bound parts, either via projection onto analytically or 
numerically known states or via spatial separation of the 
ionizing part at large distances on the grid. After removal 
of the bound parts we propagate the remaining ionizing 
part of the wavefunction backwards in time without tak- 
ing account of the interaction with the light field using 
two different Hamiltonians, once including the potential 
V(r): 



9 T 



{ion \r,t)=(\+V{r) 



(ion) 



(r,t), (9) 
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and once as a free particle: 



(10) 



In order to calculate the time delay At# 4 .r for a given 
region R, the wavcpacket has to be located outside of R 
at the start of the back-propagation and the propagation 
needs to be terminated as the wavepacket reaches the 
center of R, i.e. the location of the residual target ion. 
The latter point will be further discussed in the applica- 
tion of the method below. 



III. APPLICATION TO SINGLE 
PHOTOIONIZATION BY AN XUV PULSE 



TABLE I: Results of numerical calculations for the times 
ty g ,R, ^(o) R and the time delay At* si _H for different spatial 

steps Sx and a fixed time step of 8t = 0.002. Results are ob- 
tained for ionization from the ground state of the ID Yukawa 
potential and R = [0, ±460]. The parameters of the XUV 
pulse were: peak intensity I = 1 X 10 15 W/cm 2 , frequency 
ui = 100 eV, pulse duration r = 400 as, and carrier-envelope 
phase (f> = 0. 



A. Model systems 



spatial step 








0.5 


275.3389 


274.4569 


0.8820 


0.2 


240.6028 


241.5927 


-0.9899 


0.1 


236.8736 


237.9043 


-1.0307 


0.05 


235.9751 


237.0109 


-1.0358 


0.02 


235.7258 


236.7628 


-1.0370 


0.01 


235.6903 


236.7274 


-1.0371 



The theoretical method outlined above is, in general, 
applicable to ionization of an atom or molecule in any 
light field. Here we present results for the application 
to photoionization of an electron initially bound in two 
different model potentials. First, we used a short-range 
Yukawa potential in ID: 



V Y (x) 



Z 



\fx 2 + . 



(11) 



where Z is the effective nuclear charge, a is the soft- 
core parameter, and b is a parameter that determines the 
effective range of this ID-potential. For our simulations 
we chose Z = 3.0, a = 2.0, and b = 30.0, which relate to 
energies of —1.6742 a.u. and —1.0124 a.u. of the ground 
and first excited states. As a long-range interaction we 
made use of the Coulomb potential in ID: 



V c {x) 



(12) 



For Z = 3.0 and a = 2.0, the energies of the lowest two 
states are —1.7117 a.u. and —1.0807 a.u., which are close 
to the energies of the Yukawa potential. 

For the interaction with the XUV light pulse we used 
length gauge, i.e. 



Might (t) = Exuv(t)x 



(13) 



where -Exuv represents a linearly polarized pulse with a 
sin 2 envelope, i.e. 



-E'xuv(i) = E sin 2 (7rt/r) sm(cot - 



(14) 



where Eq is the peak amplitude, r is the pulse duration, 
lu is the central frequency, and <j) is the carrier-envelope 
phase (CEP). 

To solve the corresponding TDSE, we used the com- 
mon Crank-Nicolson method in a grid representation. In 
general, we used a spatial step of Sx = 0.02 and a time 
step of St = 0.002 in our simulations. The grid extended 



from —4000 a.u. to 4000 a.u. for the numerical simula- 
tions of the model systems interacting with an XUV pulse 
to remain the full wavefunction on the grid. The initial 
ground and first excited states were obtained by imagi- 
nary time propagation method. We continued the propa- 
gation of the wavefunction after the interaction with the 
XUV pulse until the ionizing parts of the wavepacket 
reached a distance beyond ^ 500 and hence were 
well separated from the remaining bound parts. This 
allowed us to remove the latter parts from the grid and 
remain the ionizing parts of the wavefunction only. We 
then propagated the ionizing parts at negative and posi- 
tive x backwards in time independently, either under the 
influence of the potential, Vy or Vc, or as a free parti- 
cle. We determined the corresponding times and 
t^(o) R for both parts of the ionizing wavefunction and 

added the two contributions. In the ID calculations we 
defined the region as R = [±Xi nnor , ±x out er], where Xi nncr 
and x ou t cr < 500 are the inner and outer boundaries, 
respectively, and the ±-signs apply to back-propagation 
of the two parts of the ionizing wavepacket along the 
positive/negative a;-axis, respectively. We absorbed the 
wavefunction beyond the inner boundary dinner using the 
exterior complex scaling method [ID, [26[ . 



B. Boundaries and grid parameters 

Based on the results of recent observations 0, 0] and 
calculations 0, [Tol - |T3 . |22| , we expect that the time delay 
At<s, it R as well as the difference in the time delays for the 
ionization from different initial states AT($j, ^fj, R) are 
of the orders of a few tens of attoseconds. Resolution of 
such small times requires sufficiently small time and spa- 
tial steps in the numerical simulations. In Tables U and 
HT1 we present a set of numerical results obtained for dif- 
ferent Sx and St in the case of the Yukawa potential. We 
see that a convergence of the time delay At^^R within 
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TABLE II: Results of numerical calculations for the times 
t& g ,R an d t j (o) „ and the time delay At* r for different time 
steps St and a fixed spatial step of 8x — 0.02. All the other 
parameters were the same as in Table I. 



time step 








0.1 


237.0072 


238.0431 


-1.0359 


0.05 


236.0391 


237.0757 


-1.0366 


0.02 


235.7752 


236.8121 


-1.0369 


0.01 


235.7378 


236.7747 


-1.0369 


0.005 


235.7284 


236.7654 


-1.0370 


0.002 


235.7258 


236.7628 


-1.0370 


0.001 


235.7254 


236.7624 


-1.0370 



less than 0.001 a.u. (i.e. < 0.025 as) is reached for a time 
step of St = 0.002 and a spatial step of Sx = 0.02. Sim- 
ilar conclusions hold for our studies with the Coulomb 
potential as well. 

As mentioned above, the time delay At^ it R depends 
on the size of the region R, and should be negative and 
converge to a finite limit for short-range potentials only. 
To test these expectations, we performed a set of simu- 
lations for the time delays for photoionization from the 
ground and excited states of both potentials as a func- 
tion of the outer boundary x uter by fixing Xinncr = 0, i.e. 
for R = [0,x outC r]- As expected, the values for the time 
delays are negative and decrease for an increase of the 
region R for each of the results presented in Fig. [T] For 
the Yukawa potential (panel a) convergence is found for 
outer boundaries Xoutcr > 150. Consequently, for large 
values of the outer boundary we obtain a well-defined 
value for the time difference AT of the time delays for 
ionization from the ground and the excited states. 

In contrast, our results do not show a convergence for 
the time delays as a function of the outer boundary in 
case of the long-range Coulomb potential (see Fig. D^b)). 
This reflects the well-known logarithmic divergence of the 
time delay for this kind of potential and, hence, for ion- 
ization from any bound state within the potential. Of 
course, in these cases a WS time delay as the derivative 
of the phase shift (c.f. Eq. ((4])) cannot be defined as well, 
since its derivation requires a finite limit of Ai 1 j, i ./j_ i . 0O . 
ft is interesting to point out that the results in Fig. [ljb) 
further show that the logarithmic divergence is, in gen- 
eral, still present for the difference between a pair of time 
delays obtained for the ionization from two different ini- 
tial states. Thus, such a time delay difference docs not 
simply depend on the short-range character of the poten- 
tial but contains information about the long-range part 
of the potential and is therefore not well-defined as well. 
The present results agree well with the conclusions of 
early works on time delays (24[. However, our results for 
the time delays for any finite region are finite and there- 
fore allow us to study certain aspects with respect to the 
parameters of the XUV pulse and the effects of an IR 
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FIG. 1: (Color online) Time delays Aty it R and time difference 
AT($j, tyj, R) as a function of the outer integration bound- 
ary x ou tcr for two potentials: (a) short-range Yukawa poten- 
tial and (b) long-range Coulomb potential. Time delays ob- 
tained for the ground and first-excited states are represented 
by blue dashed lines and green dash-dotted lines, respectively; 
while the red dotted lines show the results for the time dif- 
ference between the delays. In (b) the black solid and blue 
dashed lines correspond to two different forward propagation 
distances: (xf orward ) = 2000 and 3000, respectively, for the 
ionization from the ground state. In all calculations we have 
used an XUV pulse with peak intensity I = 1 x 10 15 W/cm 2 , 
central frequency uj — 100 eV, pulse duration r = 400 as, and 
carrier-envelope phase <f> = for the ionization. 



streaking in the Coulomb case. 

It is necessary to point out that for the Coulomb po- 
tential the time delay also depends on the distance that 
the ionizing wave packet is propagated in forward di- 
rection. This is due to the long-range character of the 
Coulomb potential since the central momentum of the 
ionizing wave packet decreases with increase of the for- 
ward propagation distance. Thus, the velocity of the free 
particle during the back propagation decreases as well. In 
Fig. QJb) we show this effect by presenting results for the 
time delay from the ground state for two forward propa- 
gation distances: (xf orW ard) = 2000 (black solid line) and 
3000 (blue dashed line). As expected, the time delays 
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FIG. 2: (Color online) Time delays and difference between 
time delays as a function of inner integration boundary x'i nncr . 
Symbols and laser parameters are the same as in Fig. [T] We 
also plotted the WS time delays as black dots in this figure. 



method up to \x\ — 500, projected the numerical solution 
onto the appropriate plane-wave solutions for x — > ±oo, 
and obtained the WS time delay for the scattering pro- 



cess At 



(scat) 
WS 



as the derivative of the phase shift in the 



plane wave propagating in positive x-direction with re- 
spect to the energy of the incident particle. In order to 
take account of the energy spread of the ionizing wave 
packet in a specific photoionization process, we averaged 
Atfc? over the energy spectrum of the wave packet, 
as obtained in our time-dependent numerical simula- 
tions. Finally, we considered the photoionization as a 
half-scattering process and divided the result of the av- 
erage by two. The resulting WS time delays for pho- 
toionization are shown as black dots in Fig. [2] and are 
in good agreement with our numerical results, obtained 
from the time-dependent calculations, for Xi nnor = and 
^outcr = 500. This is in support of the applicability of our 
approach to obtain time delays from the time-dependent 
numerical simulations. 



for (xforward) = 3000 are slightly smaller than those for 
(^forward) = 2000. However, this small dependence on the 
forward propagation distance does not change our con- 
clusions regarding the convergence of the results towards 
infinite regions. 

We also note from the results in Fig. Q] that the time 
delay increases most strongly in the region close to the 
center of the potential, where the potential changes dra- 
matically. This indicates that the results should depend 
on the choice of the inner boundary x; nnor of the region 
R. To study this feature, we fixed the outer boundary of 
R at a; ou tor = 500, which is large enough to obtain con- 
verged results in case of the Yukawa potential, and then 
varied the inner boundary dinner- The results in Fig. [5] 
show the expected dependence on the choice of Xi nner , 
the absolute values of the time delays decrease by half as 
dinner increases from a.u. to 5 a.u.. In the remainder of 
the present studies we have chosen the Xj nncr = a.u. as 
the inner boundary, since this value corresponds to the 
expectation value of x for all the bound states investi- 
gated here. 

C. Wigner-Smith time delay for photoionization in 
a short-range potential 

In case of the Yukawa potential, wc also compared our 
numerical results, obtained in the time-dependent numer- 
ical simulations, with calculations of the WS time delay 
as a derivative of the induced phase shift, c.f. Eq. (|4|), 
obtained from a time-independent scattering approach. 
In order to obtain the latter for photoionization by a 
light pulse with finite duration, we first considered the 
scattering of an electron, incident from x = — oo with a 
momentum k, of the short-range Yukawa potential. We 
solved the corresponding time-independent Schrodingcr 
equation numerically using the fourth order Runge-Kutta 



D. Dependence of time delay on XUV pulse 
parameters 

Next, we studied the dependence of the time delay on 
the parameters of the XUV ionizing pulse for photoion- 
ization from the ground state of the Yukawa potential. In 
Fig. [3]we present our results as functions of (a) the XUV 
frequency at a fixed pulse duration of r = 400 as and (b) 
the duration of the XUV pulse at a fixed frequency of 
u = 100 cV. The peak intensity was I = 1 X 10 15 W/cm 2 
and the CEP was <f> = in each of these simulations. 

The results agree well with qualitative expectations. 
The absolute value of the time delay decreases towards 
zero as the frequency of the ionizing XUV pulse and, 
hence, the final kinetic energy of the emitted electron in- 
creases (cf. panel a). This is due to the fact that the effect 
of the potential on the motion of the electron becomes 
negligible in the limit of infinite large kinetic energy of 
the electron (i.e. infinite large XUV frequency) and there- 
fore the time spent in the potential approaches that of 
the free particle in this limit. 

We further find that the absolute value of the time de- 
lay decreases with an increase of the XUV pulse duration 
(panel b). This dependence is closely related to that pre- 
sented in Fig. GDJa) and can be qualitatively understood 
as follows. Due to the finite pulse duration the ionized 
electron wave packet has a certain bandwidth about a 
central kinetic energy. Consequently, the time delay ob- 
tained for the wave packet can be considered as an av- 
erage over contributions at particular electron energies 
within the bandwidth (weighted by the ionization prob- 
ability at a given energy). As indicated by the results in 
panel (a) the time delay does not change linearly with 
the kinetic energy. Therefore, the time delay obtained 
for a wave packet will be smaller than its contribution 
at the central kinetic energy or the expectation value of 
the kinetic energy. This difference decreases and, thus, 
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FIG. 3: (Color online) Time delays for ionization from the 
ground state of the Yukawa potential as functions of (a) the 
XUV photon frequency (r = 400 as) and (b) the pulse dura- 
tion of the XUV pulse (w = 100 eV). Comparison between re- 
sults from the present TDSE calculations (red diamonds) and 
those for the Wigner Smith time delay (blue open squares) is 
shown. Other laser parameters were: / = 1 x 10 15 W/cm 
and 6 = 0. 



the time delay for the wave packet increases as the en- 
ergy bandwidth of the wave packet decreases, i.e. as the 
pulse duration increases. Furthermore, it is found that 
the expectation value of the kinetic energy of the ion- 
izing wavepacket increases as the XUV pulse duration 
increases, which also causes the time delay to increase. 

For each of the results from the time-dependent sim- 
ulations presented in Fig. [3] we also calculated the WS 
time delay for photoionization, as described in the pre- 
vious subsection. Again, we found excellent agreement 
between the results from the time-independent scatter- 
ing approach (blue open squares in Fig. [3]) and our nu- 
merical time-dependent simulations. We may note par- 
enthetically, that our simulations results also agree well 
with those of classical calculations (not shown) in which 
the time delays are determined by 



At 



classical 



i 



-dx- 



^outcr 



o ^2((E) - V(x)) y/2(E) 



(15) 



where (E) is the expectation value of the kinetic energy 
E of the electron for a given ionized wave packet. 



IV. STREAKING OF PHOTOIONIZATION 
PROCESSES BY NEAR-INFRARED FIELD 

In an attosecond streaking experiment [27[ , an IR field 
is used to map time information to the momentum space. 
One of the questions that arise in this context is whether 
or not the streaking field influences the observed quan- 
tities. In this section we therefore apply our theoretical 
method to such a situation. We will first discuss how 
the streaking field can be included in the numerical sim- 
ulations of time delays and then study the impact of a 
streaking IR field on the time delays for the short-range 
Yukawa as well as the long-range Coulomb potential by 
varying the parameters of the streaking field. 



A. Streaking field in the numerical simulation 

The streaking field is represented by an additional po- 
tential V^trcak which we consider as part of the poten- 
tial V(r) in Eq. ©. Thus, in the present calculations 
the streaking field is considered on equal footing with 
the atomic potential. After (forward) propagation of the 
wavefunction from its initial state and separation of the 
bound and ionizing part of the wavefunction, we then 
propagate the ionizing part of the wavefunction back- 
wards once within the combination of the atomic poten- 
tial and the streaking field and once as a free particle. 

As a result we obtain the time delay, associated with 
the ionizing part of the wavefunction in the combined 
potential of the short- or long-range interaction and the 
streaking field as 



At 



(IR) 



t 



(IR) 



*(0) B 



(16) 



where , R is the time the ionizing wave packet spends 
in region R in the presence of the (short- or long range) 
potential and the IR streaking field, and t^to) R is the 
time for the free particle case. 

Here, the forward propagation of the wavefunction 
has to be continued as long as the IR streaking field is 
present. As noted above, for the long-range Coulomb 
potential the time delay depends on the distance the 
wavepacket is propagated in forward direction. In or- 
der to keep the corresponding error small in our current 
analysis we used a large grid of —13000 a.u. to 13000 a.u. 
and stopped the forward propagation when the expecta- 
tion value of the ionizing wave packet reaches 8000 a.u.. 
We increased the spatial step to Sx = 0.1 and the time 
step to 8t = 0.02 as compared to the previous calcula- 
tions. Test calculations showed that the relative error of 
the present results is about 1%. In order to compare with 
the results, presented above, we considered the same pa- 
rameters for the Yukawa and the Coulomb potential as 
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FIG. 4: (Color online) Time delays (upper row) and time de- 
lay differences (lower row) as a function of outer boundary 
Pouter of R for Yukawa potential (left column) and Coulomb 
potential (right column) . For each potential we have centered 
the XUV pulse at two different positions which correspond to 
the maximum (blue dash-dotted line) and zero (green solid 
line) of the IR vector potential respectively. The XUV pa- 
rameters are: Jxuv = 1 X 10 15 W/cm 2 , wxuv = 100 eV, 
Txuv = 400 as, and 4>xuv = 0. The IR parameters are: 
7ir = 1 X 10 12 W/cm 2 , Air = 800 nm, N m = 3 cycle, and 
(f>m. = 0. The small box in (d) shows the long-range behavior 
of the two curves. 



before. We chose the first excited state as the initial state 
and checked that the ionization induced by the IR field 
is negligible up to an streaking field intensity of 1 x 10 13 
W/cm 2 , which is large enough for the streaking purpose. 



B. Effects of the probing pulse on time delays 

In order to analyze the effect of the IR probing field 
we obtained the time delay for photoionization from the 
first excited state of the Yukawa as well as the Coulomb 
potential in the streaking field. To this end, we applied 
the ionizing XUV pulse centered at the maximum of the 
IR streaking field (zero of the vector potential) as well 
as centered at the central zero of the IR streaking field 
(maximum of the vector potential) . In the upper row of 
Fig- SI the results for the time delays in the streaking field 
(green solid lines: XUV centered at zero of the vector 
potential; blue dashed lines: XUV centered at maximum 
of vector potential) are shown as a function of the outer 
boundary x ou t e r of the region R for the Yukawa (left) and 
the Coulomb potential (right). As in the results without 
streaking field, we see that there is a well-defined limit 
of the time delays for the short-ranged Yukawa potential 
as the region R increases, while there is no convergence 
found for the Coulomb potential. 

In order to see the effect of the probing field, we present 



in the lower row of Fig. [4] the difference between the time 
delays in the streaking field to those without streaking 
field as a function of the outer boundary of the region R, 
i.e. 



AT = At 



(IR) 



At 



(17) 



Although there is no well-defined limit of the time delays 
for infinite regions in the Coulomb case, neither with nor 
without streaking field, for any finite region the time de- 
lays are well-defined and the effect of the streaking field 
can be analyzed. The same argument applies to the weak 
dependence of the Coulomb results on the distance of for- 
ward propagation in our simulations. 

For both potentials, we see that the time delay dif- 
ference AT oscillates for neuter < 400. This oscillation 
is due to the presence of the IR field, since the ionized 
wave packet propagated up to about x ~ 400 before the 
IR streaking field ceased in the present simulations. We 
note that the differences AT arc small, less than 3% for 
the Yukawa potential and less than the numerical error 
of 1% for the Coulomb potential, compared to the time 
delays induced by the atomic potentials themselves. In 
the present calculations we therefore do not find a signifi- 
cant effect of the streaking field, neither for a short-range 
nor for a long-range potential. This is in contrast with 
the findings of some previous theoretical work [t| [H], [22| . 

Before we continue to further study the influence of the 
IR streaking field, we note a subtle point in the results 
obtained for the Coulomb potential, which are presented 
in Fig. 0] While neither the time delays with and with- 
out streaking field converge as a function of the outer 
boundary x utcr, we find a converged result (within the 
numerical error) for the time delay difference AT, if the 
XUV pulse is centered about the zero of the vector po- 
tential of the streaking field (sec green solid line in Fig. 
Etd)). This occurs since in this case the momentum dis- 
tribution of the ionizing wave packet at the end of the 
forward propagation is the same as that of the no streak- 
ing field case. In contrast, if the XUV pulse is applied 
at the maximum of the vector potential of the streak- 
ing field, the final momentum distribution is shifted and 
thus no convergence of the time delay difference within 
the range of present boundaries is found (see blue dashed 
line in Fig.gjd)). 

The conclusion that the streaking field does not influ- 
ence the time delay significantly holds over a large range 
of XUV frequencies as well as for intensities of the streak- 
ing field up to about 10 13 W/cm 2 . In Fig. Owe present 
the results for the time delay obtained in (a) the Yukawa 
potential and (b) the Coulomb potential with (blue open 
squares) and without (blue filled circles) streaking field 
as a function of the XUV frequency. Since the results 
are in close agreement we also show the relative differ- 
ence between them (green filled triangles), which does 
not exceed 5% and 2% in the Yukawa and Coulomb case, 
respectively. 

As one would expect, the relative difference between 
the results for the time delay obtained with and with- 
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FIG. 5: (Color online) Time delays with (blue open squares) 
and without (blue filled circles) IR streaking field as well as 
relative differences between the results (green filled triangles) 
as a function of XUV central frequency for (a) Yukawa poten- 
tial and (b) Coulomb potential. The XUV pulse was centered 
at the middle (maximum vector potential) of the IR field. 
Other laser parameters are the same as in Fig. [4] For the 
Coulomb case the time delays are calculated at x ou tcr = 800. 



out streaking field do increase with an increase of the IR 
streaking field intensity. This can be clearly see from the 
results shown in Fig. [6j It appears that for IR inten- 
sities up to 10 12 W/cm 2 the relative difference between 
the results is small enough such that there is no signif- 
icant effect on the time delay. While the relative dif- 
ference quickly increases beyond 10% in the case of the 
Yukawa potential with a further increase of the IR inten- 
sity, the 10%-limit is not reached for an IR intensity of 
10 13 W/cm 2 in case of the Coulomb potential. 



V. CONCLUSION 



0.3 




IR intensity (W/cm 2 ) 



0.08 




IR intensity (W/cm 2 ) 



FIG. 6: (Color online) Relative differences of time delay as 
function of IR intensity for (a) Yukawa potential and (b) 
Coulomb potential. The XUV pulse is centered at the middle 
of the IR field. Laser parameters are the same as in Fig. [4] 
expect /ir is changing. For the Coulomb case the time delays 
are taken at x ou tcr = 400, 800 and 1200. 

tive of the phase shift with respect to the energy, for the 
short-range potential. In contrast, the numerical results 
in case of the Coulomb potential are finite for any finite 
region, but they do not converge as the region increases 
to infinity, as expected. The well-defined time delays for 
a finite region enabled us to study the impact of a near- 
infrared streaking (or probing) pulse for both potentials. 
Our results show that the effect is small as long as the 
intensity of the probing field is below 10 13 W/cm 2 . 



In summary, we have applied a fundamental definition 
of time delay to time-dependent numerical simulations 
on the grid. To this end, we have obtained the difference 
between the time a particle spends in a finite region of a 
potential and the time a free particle spends in the same 
region using a back-propagation technique. The method 
is applied to photoionization of an electron in a short- 
range Yukawa as well as a long-range Coulomb potential 
by an attosecond XUV pulse. It is found that the nu- 
merical results are in excellent agreement with those for 
the (asymptotic) WS time delay, obtained as the deriva- 
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